Azimuthal anisotropy of direct photons 
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The electromagnetic bremsstrahlung produced by a quark interacting with nucleons or nuclei is 
azimuthally asymmetric. In the light-cone dipole approach this effect is related to the orientation de- 
pendent dipole cross section. Such a radiation anisotropy is expected to contribute to the azimuthal 
asymmetry of direct photons in pA and AA collisions, as well as in DIS and in the production of 
dileptons. 
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I. INTRODUCTION 

Direct photons, i.e. photons not from hadronic decay, 
are of particular interest, since they do not participate 
in the strong interaction and therefore carry undisturbed 
information about the dynamics of the primary hard col- 
lision. 

Here we present the basic color-dipole formalism for 
calculating the azimuthal distribution of direct photons 
radiated by a quark interacting either with a nucleon 
or nuclear targets. For this purpose we further develop 
the dipole approach proposed in [l|, l2|] f° r electromagnetic 
bremsstrahlung by a quark interacting with nucleons and 
nuclei. This technique can be applied to dilepton [!, 0, [f| 
and prompt photon 0, @| production in pp, pA and heavy 
ion collisions. It can be also used for calculating the 
azimuthal angle dependence in the radiation of dileptons, 
or in deep inelastic scattering. 

An azimuthal asymmetry appears due to dependence 
of the interaction of a dipole on its orientation. Indeed, 
a colorless qq dipole is able to interact only due to the 
difference between the impact parameters of q and q rel- 
ative to the scattering center. If b is the impact param- 
eter of the center of gravity of the dipole, and r is the 
transverse separation of the q and q, then the dipole in- 
teraction should vanish if r _L b, but should be maximal 
if r || b. One can see this on a simple example of dipole 
interacting with a quark in Born approximation. The 
partial elastic amplitude reads, 



lmf$ q (b, f) = 



d 2 qd 2 q'a s (q 2 )a s (q' 2 ) 
(q 2 + fi 2 )(q' 2 + ^ 2 ) 



,iq--(b+f/2) 



,iq-'-(b+r/2) _ iq'-(b-f/2) 



(1) 



Here we assume for the sake of simplicity that q and q 
have equal longitudinal momenta, i.e. they are equally 
distant from the dipole center of gravity. The general case 
of unequal sharing of the dipole momentum is considered 
later in (|23|) . We introduced in |T]) an effective gluon 



mass [l which takes into account confinement and other 
possible nonperturbative effects. 

Integrating in JTJ) with a fixed a s we arrive at, 



lm fq q (b,r) 



9 



-K [n 



b- 



(2) 

where Kq(x) is the modified Bessel function. This ex- 
pression explicitly exposes a correlation between r and b: 
the two terms cancel each other if b ■ r = 0. 



II. DIRECT PHOTONS: DIPOLE 
REPRESENTATION 



The radiation of direct photons, which in the parton 
model looks like a Compton process gq — > jq, in the 
target rest frame should be treated as electromagnetic 
bremsstrahlung by a quark interacting with the target. 
In the light-cone dipole approach the transverse momen- 
tum distribution of photon bremsstrahlung by a quark 
propagating interacting with a target t (nucleon, t = N , 
or nucleus, t = A) at impact parameter 6, can be written 
in the factorized form [2] , 



dai T ^ x {b, Pl a) 



— Y 



d 2 r 1 d 2 r 2 e 



d{lna) d 2 p d 2 b 
x <p* q {a,r 1 )<p lq (a,r 2 )F t (b,ar 1 ,ar 2 ,x). (3) 



Here 



p and a = p^/p^ are the transverse and 
fractional light-cone momenta of the radiated photon, 
4> lq {a,r) is the light-cone distribution amplitude for the 
<77 Fock component with transverse separation r, and 
F t (b, af\, ar 2 , x) is an effective partial amplitude to be 
discussed below. The product of the distribution ampli- 
tudes, summed in (J3j) over initial and final polarizations 
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of the quark and photon, reads Q, 



differential cross section, which can be calculated as, 



(j> iq {a,r 1 )4> lq {a,r 2 ) = -^rn q a 

in,f 



x ^a 2 K Q (am, q r 1 )K Q (am q r 2 ) 



+ [1 + {l-a) 2 } r —K l {am q r l )K l {arri q r 2 )X. (4) 

T\T 2 ) 

Here m q is the effective quark mass, which is in fact an 
infra-red cutoff parameter, and can be adjusted to pho- 
toproduction data Q , or shadowing , and whose value 
is m q « 0.2 GeV. 

In equation §3§ the effective partial amplitude 

Ft (b, afi , af 2 , x) is a linear combination of qq dipole par- 
tial amplitudes at impact parameter b, 



F t (b, afi,ar 2l x) 



Im 



fUb, an,x) + f q Jb, ar 2 ,x) 



f Jb,a{ri -r 2 ),x) 



(5) 



where x is Bjorkcn variable of the target gluons. 



III. AZIMUTHAL ASYMMETRY IN 
QUARK-NUCLEON COLLISIONS 

In the case of a nucleon target (t = N), the partial 
elastic amplitude fq q {b,r) of interaction of the qq dipole 

with a proton at impact parameter b, is related to the 
dipole cross section as, 



d 2 blmf q N Jb,r) 



(6) 



where cr qq (r) is the total cross section of a qq - proton 
collision. Here and further on, unless specified otherwise, 
the dipole cross section and partial amplitudes implicitly 
depend on the Bjorken variable x of the target gluons. 

The cross section <Jq q {r) has been rather well deter- 
mined by data on deep-inelastic scattering [13] • With 
this input, and using Eq. (|3|), one can calculate the in- 
clusive differential cross section of direct photon emission. 
This was done in Q for pp collisions, with results in good 
agreement with data. 

Using the partial elastic amplitude fq q (b,r) one can 
also calculate the differential elastic cross section of 
dipole- nucleon scattering. Neglecting the real part, the 
amplitude reads, 



da 



(qq)N 



(0 



dq\ 



1 

4tt 



K(r)] 2 



d 2 be^hmf q N q (b,f) 



I6tt 



exp 



-B 



(qq)N, 



(r)q 2 T ■ (7) 



In the second line of this equation we rely on the small- 
qr approximation. This defines the forward slope of the 



B 



(qq)N 



d 2 ss 2 \Tnf q N q {s,f). (8) 



The slope for small-dipole-proton elastic scattering was 
measured in diffractive electroproduction of p-mesons at 
high Q 2 at HERA [HI- The measured slope, B ( j q)N ' (r) w 

5 GeV -2 , agrees with the expected value B^ q ' N (r) ~ 
B%/2. 

The objective of this paper is the azimuthal asymme- 
try of photon radiation. First of all, we calculate the 
asymmetry of the cross section Eq. (J3J) for quark-nucleon 
collisions. The only vector available for such asymmetry 
is the impact parameter 6, and therefore we should trace 
a correlation between the vectors p and b. The popular 
correlation function is defined as, 



v?, N (b,p,a) = (u 2 )0 p 




r 2tt , , » do qT ^~' x (b, P ,a) 
Jq a<PpV2 d(lna)d 2 pcPb 
f 27r daiT^-,x {biPia — 
JO al fP d(lna)d 2 pd 2 b 



(9) 



where the averaging is performed integrating in (|3|) over 
the azimuthal angle <p p of the transverse momentum p. 



IV. RADIATION PRODUCED BY A QUARK 
PROPAGATING THROUGH A NUCLEUS 

In this case the partial amplitude to be used in |(SJ), for 
a qq dipole colliding with a nucleus at impact parameter 
b, reads, 



lmf q A q (b,r) = 1- 



l-^ q (r)f A (b,f) 



1 — exp 



- 5 <(r)f A (6,f) 



(10) 



The effective nuclear thickness Ta is defined as [I2| , 
2 



T A (b,r) 



d s Imfq Js, r) T A (b + s) . (II) 



where the nuclear thickness function is defined as an in- 
tegral of the nuclear density along the particle trajectory, 
T A{b) = J^oodzpAib^z). 

Calculating v\ (b,p, a), we can average over 



/(£±Y\ 7^ (tlY d<J qA ^ x (b,p,a) 
\\pb j ] ° C J P \pb ) d(lna)d 2 pd 2 b ' 

(12) 

analytically. Instead of integration over direction of p at 
fixed b, one can integrate over direction of b at fixed p. 
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The advantage of such a replacement is obvious: all the 
^-dependence in ([3]) is located in the effective amplitude 
F A and it has an explicit and simple form. 

Indeed, the mean value of s 2 is according to © (s 2 ) w 
0.4 fm 2 , which is much smaller than the heavy nucleus 
radius squared, R 2 A . Therefore we can expand T A (b + s) 
as 



T A (b + s) = T A (b) + ^T A (b) + \ l^j T'i(b) + ... 

(13) 

Correspondingly, the partial amplitude (|10p can be ex- 
panded as, 



bnffsfar) ~ 1 - exp 



x l-^T» 7l (6,r3 



-l a f q (r)T A (b) 



As far as the partial amplitudes (s, r) and their 
asymmetric part Eq. (|16[) are known, one can calculate 
the azimuthal asymmetry of photons radiated in quark- 
nucleus collisions, 

qA ^ da^ x (b,p,a) ^ f 2 2 



d(lna) d 2 p cPb 



inj 



x e^^- ?2 )^ g («, ri)^(a, r 2 )F^(6, ar 2 ). (19) 
where 



^ A (6,afi,af 2 ) = Im/^(6,Q!ri) + Im/^(6,Q!r2) 



(20) 



Notice that the cross section in the left hand side 
of Eq. (fP9"| can also be calculated without using the 
expansion (| 13[) . relying on the eikonal approximation, 
T A (b,r) « T A (b), which is known to be quite accurate 
for heavy nuclei. 



T» 72 (M-TA(fo) 7l 2 (M 



(14) 



V. PARTIAL DIPOLE AMPLITUDE fg (b, r) 



where 



7„(M= / d"slmm(s,f)(s-br 



(15) 



Integrating the amplitude (|T4j) together with v 2 over 
</>& we find that the first two terms in the curly brackets 
in (Ti"4")) give zero, and the rest is, 



1 

2^ 



hlmf A (b,r)v 2 (<Pb) 



T' A \b)g(r)~T A \b)h(r) . (16) 



Here 



(p-s) 2 



(17) 



h(r) = J d 2 s 1 d 2 s 2 lu,fg{s 1 ,r)lu,fg{s 2 ,T) 
(si ■p)(s 2 -p) 



(s*i • s 2 ) 



(18) 



Eq. (|16|) shows that the azimuthal asymmetry is 
strongly enhanced at the nuclear periphery. Indeed, at 
small impact parameters the amplitude Eq. (|16[) is sup- 
pressed by the factor exp[— ^a qq {r)T A {b)\, and moreover, 

T' A {b) « -2p b/^R^¥ and T» » 2p R 2 A /(R 2 A - 
6 2 ) 3 / 2 are vanishingly small and peak at the periphery 
(po ~ 0.16 fm~ 3 is the central nuclear density). The 
smallness of both the amplitude and azimuthal asym- 
metry justifies also the expansion made in (jTJ 



The next step is to model the partial dipole ampli- 
tude fq q (b, r) . An azimuthal asymmetry can only emerge 
if the amplitude f^(b,r) contains a correlation between 

the vectors b and r. If such a correlation is lacking, the 
functions Eqs. (fT7|) and (fl"8]) are equal to zero. A model 
for fq q {b, r) having no b — r correlation was proposed in 

It is rather straightforward to calculate the partial am- 
plitude within the two gluon exchange model [l4[ , 



Im/^(6,rT 



x e 



2_ f d 2 qd 2 q'a s (g 2 )a s (q' 2 ) 
Sir 2 J (q 2 

ib-(q-q') ^ _ e 'ii 

F N (q-q')-F^\q,q-') 



p 2 ){q' 2 +p 2 ) 
\ - e~* 



(21) 



where Fjy(k) = (v^aH exp(«fc • /3i)|\E r at) is the nucleon form 
factor, and F ( ^ q) (q,q') = (V N \ exp[iq ■ p x - iq' ■ p 2 ]|*iv) 
is the so called two-quark nucleon form factor. Both can 
be calculated using the three valence quark nucleon wave 
function ^ N (pi, p 2 , p 3 ). 

An effective gluon mass p is introduced in (|2ip in order 
to imitate confinement. We fix its value at p = in 
order to reproduce the large hadronic cross sections. 

The Born amplitude is unrealistic since leads to an 
energy independent dipole cross section a qq (r, x) . This 
dipole cross section has been well probed by measure- 
ments of the proton structure function at small Bjorken 
x at HERA, and was found to rise towards small x, with 
an x dependent steepness. In fact, it can be expressed 
via the unintegrated gluon density T(x, g 2 ), 



*(r,*) = 



a s (q 2 )T( Xl q 2 ). (22) 
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Analogously, the partial amplitude for dipole-nucleon 
clastic scattering at impact parameter b between the cen- 
ters of gravity of the dipole and nucleon reads, 



-iq-r(3 _ J,q-r{ 



-«) (, 



d 2 qd 2 q' 
q 2 q' 2 

iq -rf3 



a s F(x,q,q'y b ^"> 



-i<f-r(l-l3) 



(23) 



Here the dipole has transverse separation r, fractional 
light-cone momenta of the quark and antiquark, 1 — (3 
and P respectively. Since the radiated photon takes away 
fraction a of the quark momentum, the corresponding 
dipole has j3 = 1/(2 — a). The impact parameter b of 
the dipole is the transverse distance from the target to 
the dipole center of gravity, which is shifted towards the 
fastest q or q in accordance with 



In (f2"3"|) a s = y / a s (q 2 )a s (q' 2 ), and we introduced the 
off-diagonal unintegrated gluon density T{x, q, q'), which 
in the Born approximation limit takes the form, 



T(x,q,q) 



IT 



F N (q-q') 



F^\q,q') .(24) 



Besides, the partial elastic amplitude Eq. (f2"3"| . should 
satisfy the conditions Eqs. (J6]) and ([8]). For the dipole 
cross section we rely on the popular saturated shape fioj ] 
fitted to HERA data for F%{x,Q 2 ) and we choose the 
following form of T(x, q, q' ), 



F(x, q, q") 



3o-p 
16 7r 2 a s 
1 



X exp^--i?2( x ) {q 2 +q r - 
e*v[~R 2 N {q-q') 2 /2\ , 



(25) 



where er = 23.03 mb, FUx) = 0.4 fm x (x/xq) ' 144 
with x = 3.04 x 10~ 4 and i = pjjl 7]. We 
assume here that the Pomeron-proton form factor has 

so the 

•pp _ 



the Gaussian form, Fjp(k^) — exp(—k^R N /2), 
slope of the pp elastic differential cross section is B e[ 
2R%- 



2a'jp ln(s/so), where a'jp ~ 0.25 GcV " is the slope 



of the Pomeron trajectory, sq = 1 GeV . 



c \)/3 



is the part of the slope of elastic cross section related to 
the Pomeron-proton form factor. 

With this unintegrated gluon density the partial am- 
plitude Eq. ([23|) can be calculated explicitly, 



CO 



8irB e 



exp 



[b + r(l-P)} 2 
2B el 





(b~r/3) 2 ~ 




r 2 


+ exp 


2B el 


— 2 exp 


Rq 



[b+(l/2-P)r\ 2 
2B e , 



>, (26) 



where B e i(x) — R 2 N + Rq(x)/8. This amplitude satisfies 
the conditions Eqs. ([6|) and ((8]). This expression also goes 
beyond the usual assumption that the dipole cross section 
is independent of the light-cone momentum sharing (3. 
The partial amplitude Eq. (|26|) does depend on f3, but 
this dependence disappears after integration over impact 
parameter b. 



VI. NUMERICAL RESULTS 

Now we are in a position to calculate v\ (b,p, a). Ex- 
amples of quark- nucleon collisions radiating a photon, 
with a = 1 and at different impact parameters and en- 
ergies, are depicted in Fig. [TJ The results show that 
the anisotropy of the dipole interaction rises with im- 
pact parameter, reaching rather large values. As func- 
tion of the transverse momentum of the radiated pho- 
tons, t^^p, a) vanishes at large px- Such a behavior 
could be anticipated, since the interaction of vanishingly 
small dipoles responsible for large p is not sensitive to 
the dipole orientation. 

The next step is calculating the azimuthal asymmetry 
in quark- nucleus collisions. The results are plotted in 
Fig. [3 as function of transverse momentum, at different 



impact parameters and at the energies of RHIC and LHC. 

The first observation is the smallness of v\ A , which 
is suppressed an order of magnitude compared to v% N . 
At first glance this might look strange, since the quark 
interacts with nucleons anyway. However, a quark prop- 
agating through a nucleus interacts with different nucle- 
ons located at different azimuthal angles relative to the 
quark trajectory. Their contributions to v\ A tend to can- 
cel each other, restoring the azimuthal symmetry. Such 
cancellation would be exact if the nuclear profile function 
Ta(6) were constant. We have a nonzero, but small v\ A 
only due to the variation of Ta with 6, i.e. the presence 
of finite first and second derivatives, as was derived in 

Eq. mi). 

The results of a numerical integration (without expan- 
sion (TT4"]) ). depicted in Figs. [T][21 also confirm the antici- 
pation based on Eq. (|16p that the azimuthal asymmetry 
is enhanced on the nuclear periphery. 

We used the Woods-Saxon parametrization for nuclear 
density [l5| . The anisotropy of electromagnetic radia- 
tion appears only on the nuclear periphery and accord- 
ing to (Unj) is extremely sensitive to the behavior of the 
nuclear thickness function at the very edge of the nu- 
cleus. Electron scattering data, which is the main source 
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FIG. 1: The anisotropy parameter v\ (b,p,a) as function of 
p calculated at a — 1 for different impact parameters b and 
energies: ^/i = 200 GeV (solid, b = 0.2, 0.4, 0.6, f fm), and 
v/s = 5500 GeV (dashed, b = 0.6, 1 fm). 
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FIG. 2: Azimuthal anisotropy of direct photons with a = 1 
from quark-lead collisions at different impact parameters as 
is labeled in the plot. Solid and dashed curves correspond 
to the energies of RHIC {b = 5, 6, 7, 8, 9.5 fm), and LHC 
(b = 6, 8fm) respectively. 

of information about the electric charge distribution in 
nuclei, is not sensitive to the neutron distribution, which 
is known to be enlarged on the periphery. Therefore 
the details of the shape of the density distribution on 
the nuclear surface are poorly known. As a simple esti- 
mate of the theoretical uncertainty related to this prob- 
lem one can use an alternative parametrization of the nu- 
clear density, such as the simple and popular hard sphere 
form, p(r) = pqQ(Ra — r). We compare in Fig. [3] the 
anisotropy parameters u| (P> b, a) calculated with hard 



FIG. 3: Azimuthal anisotropy of direct photons with a = 1 
from quark- lead collisions at b — 5 and 7fm. Solid and 
dashed curves are calculated with Woods-Saxon (WS) and 
hard sphere (HS) parametrizations of nuclear density respec- 
tively. 

sphere (dashed) and Woods-Saxon (solid) parametriza- 
tions. As one could expect, the hard sphere density leads 
to a quite larger anisotropy, since the derivatives of the 
nuclear profile function are much sharper. 

VII. SUMMARY 

Summarizing, we extended the dipole description of 
electromagnetic radiation [l], Q in quark nucleon and nu- 
cleus collisions to calculation of the azimuthal angle dis- 
tribution. This problem involves more detailed features 
of the dipole amplitude, namely its dependence on dipole 
size and impact parameter, as well as on their correlation. 
We propose a simple model generalizing the unintegrated 
gluon density fitted to HERA data for the proton struc- 
ture function to an off-diagonal gluon distribution. The 
latter satisfies all the imposed boundary conditions. 

The developed theoretical tools can be applied to the 
calculation of the azimuthal asymmetry in DIS and in 
Drell-Yan reactions on a proton, as well as to the pro- 
duction of direct photons and Drell-Yan pairs in proton- 
nucleus and heavy ion collisions. 
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